March 27, 2013 11:18 WSPC/INSTRUCTION FILE oo-review-fmal 



International Journal of Modern Physics B 
© World Scientific Publishing Company 



Elementary Excitations due to Orbital Degrees of Freedom in Iron 

Based Superconductors 



WEI-CHENG LEE 

Department of Physics, University of Illinois, 1110 W. Green Street 
Urbana, Illinois 61801, United States of America 
leewc @illinois. edu 



WEICHENG LV 

Department of Physics and Astronomy, University of Tennessee 
Knoxville, Tennessee 37996, United States of America 
Materials Science and Technology Division, Oak Ridge National Laboratory 
Oak Ridge, Tennessee 37831, United States of America 

HAMOOD Z. ARHAM 

Department of Physics, University of Illinois, 1110 W. Green Street 
Urbana, Illinois 61801, United States of America 



One central issue under intense debate in the study of the iron based superconductors 
is the origin of the structural phase transition that changes the crystal lattice symmetry 
from tetragonal to orthorhombic. This structural phase transition, occurring universally 
in almost every family of the iron-based superconductors, breaks the lattice C4 rotational 
symmetry and results in an anisotropy in a number of physical properties. Due to the 
unique topology of the Fermi surface, both orbital- and spin-based scenarios have been 
proposed as the driving force. In this review, we focus on theories from the orbital- 
based scenario and discuss several related experiments. It is pointed out that although 
both scenarios lead to the same macroscopic phases and are not distinguishable in bulk 
measurements of the thermodynamic properties, the elementary excitations could be 
fundamentally different, and provide us with the possibility to resolve this long-standing 
debate between orbital- and spin-based theories. 

Keywords: orbital ordering; non-Fermi liquid; nematicity. 



1. Introduction 

Understanding the nature of the excitations in both the normal and the super- 
conducting states has been a heavily studied research topic in high-temperature 
superconductors, for it might hold the key to the pairing mechanism for unconven- 
tional superconductivity, the "Holy Grail" of the field. It is no doubt that the spin 
physics plays a crucial role in this aspect due to the presence of some magnetically 
ordered states in the parent compounds. After the doping is introduced, the super- 
conductivity emerges with the disappearance of the 'parent' magnetic order. Since 
this sort of phase diagram is so universal, it is reasonable to believe that the spin 
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fluctuations originating from the 'parent' magnetic order persist even as the long- 
ranged order is killed, and they might take a part in the formation of the Cooper 
pairs. A number of experiments support the ubiquity of the spin fluctuations among 
cuprates, heavy fermions, and iron pnictides P 

The new high temperature superconductor discovered in 2008, the iron pnictide, 
has received a tremendous amount of attention in the past few years!^ 1 ^ At first 
sight, this category of materials seems to be incompatible with superconductivity. 
The key element, Fe 2+ , has six electrons in its 3d orbitals, and Hund's coupling 
should favor ferromagnetism (parallel spins) instead of spin-singlet superconductiv- 
ity (anti-parallel spins). Indeed there were some theoretical proposals of the iron 
pnictide being a spin-triplet p-wave superconductor SI but this was ruled out very 
quickly by experiments. The hybridizations between the d orbitals of the Fe atom 
and the p orbitals of the As atom result in a strong suppression of the effect of 
Hund's coupling.^ Consequently, the magnetic moment of the Fe atom is surpris- 
ingly small, which prevents the occurrences of ferromagnetism as well as spin-triplet 
superconductivity. Experimentally it has been found that the Fe atom has a mag- 
netic moment ranging from 0.5 to 3.3 hb in different families of the iron pnictides, 
depending on the details of their electronic structures. This sensitivity of the Fe mo- 
ment on the electronic structures reflects the hybridization effect discussed above. 

Another crucial consequence due to the strong hybridization is the m ultiorbital 
nature. As seen in density functional theory (DFT) calculations J2138] ^he Fermi 
surface of the iron pnictides dissects into different pockets centered around high 
symmetry points of the Brillouin zone T and M. A closer look into the band struc- 
tures makes it obvious that the eigenstates on different parts of the Fermi surface 
have very different weights in the Fe d orbitals. In fact, it is just this multiorbital 
nature that makes the iron pnictides distinct from the cuprates. 

Fig. Q] shows a generic phase diagram of the iron-based superconductors. Let us 
first discuss the part that is remarkably similar to that of the cuprates. The parent 
compound has a stripe-like antiferromagnetic (AFM) order with the ordering wave 
vector Q — (ir,0), while superconductivity appears as the long-ranged AFM order 
is suppressed by either doping or pressure. However, a new ingredient that is not 
universally present in the cuprates is the emergence of a structural phase transi- 
tion, across which the crystal symmetry goes from tetragonal to orthorhombic. The 
structural phase transition spontaneously breaks the lattice C4 rotational symme- 
try and always occurs at a temperature (T$) which is the same as or higher than 
the magnetic transition temperature (Tn). This unusual inequality (Tg > T/v) has 
motivated theorists to propose cither orbital or spin fluctuations as the effective 
cause. The difficulty to resolve this controversy is that both theories end up with 
the same macroscopic states breaking exactly the same C4 symmetries, so that they 
arc indistinguishable in most bulk measurements of the thermodynamic properties. 

In this review, we focus on the orbital physics emerging from this multiorbital 
nature which is new and unique in the iron pnictides. It is shown that the elementary 
excitations of the orbital fluctuations is fundamentally different from those of the 
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Fig. 1. The generic phase diagram of the iron- based superconductors as a function of doping x 
and temperature T. We use Tg, TV, and Tq to label the structural, magnetic, and superconducting 
phase transition temperatures, respectively. The corresponding crystal and magnetic structures of 
each phase are displayed. 



spin fluctuations, which might give us a chance to resolve the nature of the quantum 
fluctuations from various types of spectroscopies. This paper is organized as follows. 
Sec. [5] introduces the ferro-orbital order of Fe d xz and d yz orbitals as the effective 
cause of the structural phase transition. In particular, we will show how this orbital 
order arises from both the strong-coupling and the weak-coupling limit. In Sec. |31 
the elementary excitations from the quantum fluctuations associated with the ferro- 
orbital order are discussed, and the experiments that might have detected these 
orbital fluctuations are summarized in Sec. 4. Finally, a concluding remark is given 
in Sec. 5. 



2. Orbital Order 

2.1. Structural phase transition and nematic order 

The ground state manifold of the magnetically ordered phase of the iron-based 
superconductors is 0(3) x Z 2 , where 0(3) refers to the arbitrary direction of the 
ordered moment in spin space, whereas Zi represents the two- fold degeneracy of 
the ordering wave vector (it, 0) and (0, ir). This magnetic order not only breaks the 



27, 2013 11:18 WSPC/INSTRUCTION FILE oo-review-final 



4 W.-C. Lee, W. Lv, and H. Z. Arham 

usual 0(3) spin rotational symmetry, but also breaks the C4 rotational symmetry 
of the tetragonal lattice. From a symmetry argument, accompanying the onset of 
the magnetic order, a structural phase transition could occur and reduce the crystal 
symmetry from tetragonal to orthorhombic. Furthermore, a more intriguing possi- 
bility is that the C4 symmetry breaking can arise even without breaking the spin 
0(3) symmetry. Indeed, experimentally the structural phase transition temperature 
Ts is higher than the Neel temperature Tjq in many iron-based superconductors. 
The phase sandwiched between Ts and TV is usually called the nematic phase, 
which proliferates in many strongly correlated systems P In its original meaning, 
the nematic order is characterized by a state that breaks a continuous rotational 
symmetry but remains invariant under other symmetry operations. However, in 
real materials, the underlying symmetry of the system becomes discrete due to 
the presence of the lattice. So the nematic order only breaks a discrete rotational 
symmetry, in particular the C4 rotational symmetry of the two-dimensional square 
lattice considered here. 

Understanding the origin of the C4 symmetry breaking associated with the 
structural phase transition is a fundamental question in the iron-based su- 
perconductors. Recently a lot of experimental efforts have been made to un- 
cover the ele ctronic a nisotropy in this nematic phase. Evidences from tr ansp ort 
measurement P^- ^ angl e-resol ved photoemission spectroscopy ( ARPES) and 
inelastic neutron scatterin g^ 4 | 15 | a jj SU gg es t that strong electronic anisotropy exist 
above Tjv (for a review, see Re f. ITHl) . and in some cases, it even persists above T5 
with or without uniaxial strains^-4 We note here that the orthorhombic distortion 
is very small, less than 1% of the original lattice. So it is difficult to ascribe the 
strong electronic anisotropy to lattice distortions. Therefore it is generally agreed 
that the symmetry breaking has an electronic origin, which is usually termed as the 
electron nematic order ff^l 

In this section, we focus on theories in which the orbital order is an effective 
cause of the electron nemat icity! 1 9 1 2 | 2 1] ^ g men ti onec j j n s ec [JJ iron-b ased s uper- 
conductors are intrinsically multiorbital systems, and DFT calculations^^!] show 
that all five Fe 3d orbitals have significant weights on the Fermi surface. Among five 
d orbitals, two of them, d xz and d yz , are special from the symmetry aspect. Each 
of them has only the C2 symmetry, but the system retains the C4 symmetry if they 
are degenerate. Moreover, they make dominant contributions on the Fermi surface 
compared to the other three orbitals. Therefore the minimal model for iron based 

IOQI 

superconductors is a two-orbital model with the d xz and d yz orbitals!^ 21 In such a 
model, we can define the ferro-orbital order parameter as the difference between 
the occupation numbers of the two orbitals, m = n xz — n yz . A non-zero m lifts 
the orbital degeneracy and changes sign under the C4 rotation. Therefore it can be 
identified as the electron nematic order parameter, which gives rise to the structural 
phase transition. 

What is the microscopic mechanism driving the orbital order? Before answering 
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this question, we need to address the important issue regarding the strength of 
electron correlations. Because various experiments ha ve char acterized the iron-based 
superconductors as intermediately correlated metals J24125I both strong-coupling^ 
and weak-coupling^^ theories have their own merits. So we will use arguments from 
both aspects to show how the orbital order arises in Sec. 12.21 and 12.31 respectively. 

Here we briefly mention the spin-based scenario for the structural phase transi- 
tion. In this line of thinking, the nematic order is induced by magnetic fluctuations, 
termed as the spin nematic order. Thes e theo ries can be formulated either from a 
strong-coupling J1-J2 Heisenberg or from a weak-coupling model con- 

sisting of electron and hole pockets.^3 Although starting from the opposite limits, 
both approaches end up with the same effective Landau- Ginzberg model, in which 
the nematic order can be identified as a Z2 order related to magnetic fluctuations 
centered at (7r,0) and (0,7r). Once the nematic order sets in, the magnetic fluctu- 
ations will be enhanced at a preferable wavevector without the onset of the long- 
ranged magnetic order. This approach naturally explains the proximity of structural 
and magnetic transitions. From a symmetry point of view, both the orbital order 
and spin nematic order break the same C4 rotational symmetry. Thus they are al- 
lowed to couple to each other directly. Namely, the onset of one order parameter will 
induce the other. It is therefore very difficult to distinguish their roles in causing 
the electron nematic order. However, a key point to be emphasized is that orbital 
order gives a new energy scale that is distinct from spin physics. As we will show in 
Sec. [3] and 21 orbital order leads to a new type of elementary excitations that has 
strong support from experiments. 



2.2. Orbital order from strong coupling 

The idea of orbital ord er in the insulating limit is pioneered by the seminal work 
of Kugel and Khomskii!^ We will apply their approach to the iron-based super- 
conductors, in particular a two-orbital model consisting of d xz and d yz orbitals at 
quarter filling. On each site there is one electron occupying any state that is a linear 
combination of d xz and d yz . Formally, such a state is written as 

(i) 



c( >s l —J \xz) + sm ^ 

Due to the two-fold degeneracy, we can introduce an orbital pseudospin operator T 
with T = 1/2 to represent the two orbital states, where 

lA /0 




VZ) = ( I 1 ■ (2) 



According to Kugel and KhomskiifSH this orbital degeneracy can be lifted by 
three types of interactions in a lattice system. The first one is the quadrupole- 
quadrupole interaction arising from the quadrupole moment of the 3d orbitals. In 
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terms of orbital pseudospins T, it takes the form of 

it.,. = E ■'>.:>■ >i + + { f i f ! + '■'":) ■ ( 3 ) 

ij 

Depending on the values of the exchange constants J, different orbitally ordered 
states can be stabilized. If J[j is negative and the leading energy scale of the system, 
ferro-orbital order will develop at low temperature, where m = (T?) ^ 0. 
The second one is the electron-phonon coupling, 

H ep = £ " k a\a k - E 0* + 4) Tt + ^ - | E ^ ( 4 ) 

The fist term represents free phonons with frequency ajfc, while the second term de- 
scribes the electron-phonon interaction, with gk being the coupling strength. We 
have separated the contributions from homogeneous distortions in the last two 
terms. Here 5 is the orthorhombic distortion; C is the elastic constant; and go 
is the coupling strength. We can formally integrate out the lattice variables and 
derive an effective model of pseudospins, 

Hfaty =^W? (5) 

where the exchange constants 

j.. = _^_J2^-e ik ^- r ^. (6) 

We note that only the orthorhombic distortion that couples to T z is considered 
here. As a result, the effective Hamiltonian of Eq. ([5]) takes the form of an Ising 
model. In reality, there is another phonon mode that couples to T x . Inclusion of this 
coupling will produce terms of the form 7f Tf and Tf T?, resulting in an effective 
Hamiltonian similar to the form of Eq. Nevertheless, since ferro-orbital order 
and the accompanying orthorhombic distortion are the leading instability in our 
system, it is sufficient to focus only on T z . Similar to the quadrupole-quadrupole 
interaction, electron-phonon interaction also leads to a ferro-orbitally ordered phase 
at low temperature. 

The last type of interaction is the exchange interaction. Unlike the first two 
interactions, which can be expressed exclusively by orbital pseudospins T, the ex- 
change interaction also involves the spin degrees of freedom. In many systems, it 
is the dominant energy scale for orbital and spin orders. Formally this spin-orbital 
superexchange model can be derived from a second-order strong-coupling expan- 
sion of a multiorbital Hubbard model. Considering the virtual hopping pr oces s 
d^d 1 —> d 2 d° — »• d 1 d 1 , we write down its general form along a given bond (z, j ) J32] 

^so = E E 4:i], Si , Sj A^m^ x ^w^o. ( 7 ) 

Ti,Tj Si .Sj 
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(a) f, (b) (c) (d) 




AE=-At, 2 /U AE=-4t 2 2 /U AE = -4(t, 2 +t 2 2 )/(V-J) AE = -A{t, 2 +t 2 2 )/V 

Fig. 2. Illustration of the spin-orbital coupling. We consider virtual hopping processes along the 
x direction. The intra-orbital hopping amplitudes of d xz and d yz are ti and <2, respectively. The 
energy of each orbital and spin configuration can be calculated using second-order perturbation 
theory. U, V, and J are intra- and inter-orbital Coulomb repulsion, and the Hund's coupling, 
respectively. 

where S and T are the S =1/2 spin and the T = 1/2 orbital pseudospin operators, 
respectively. Here we use (r^, t,-) and (s^ Sj) to label the orbital and spin quantum 
numbers of the intermediate high-energy state. 

Instead of presenting a formal derivation of the spin-orbital model, here we use 
simple arguments from GoodenouglPSl anc j KanamorfSl to illustrate how spin and 
orbital orders are coupled to each other. As an example, let us consider the virtual 
hopping processes between nearest neighbors (NN) along the x direction for a two- 
orbital model. We define the intra-orbital hopping amplitudes of the d xz and d yz 
orbitals as t\ and t%, respectively. The inter-orbital hopping vanishes from symmetry 
arguments. We need to consider three distinct orbital configurations as illustrated 
in Fig. [5] If both sites are occupied by the same orbitals [Fig. [^a) and (b)], the 
exchange interaction is strongly antiferromagnetic (AFM). We can write down an 
effective spin-only Heisenberg model Ti{Si, Sj) = J Si ■ Sj, where J = J\ a = 4t\/U 
if d xz occupies both sites, while J = Jib = ^k\jXJ when we have d yz orbitals on both 
sites. Here U is the intra-orbital Coulomb repulsion. In contrast, if different orbitals 
occupy the two sites [Fig.[2jc) and (d)], the exchange becomes weakly ferromagnetic 
(FM), with 

U{S i ,S i )=[ U y 2) - y^f- JSi-Sj*- U y2 Si-Si, (8) 

in the limit that the Hund coupling J is much smaller than the inter-orbital Coulomb 
repulsion V. Consequently, different orbital configurations correspond to different 
spin configurations, and vice versa. Namely, the orbital and spin orders are strongly 
coupled. If we further assume that t\ 3> <2 according to the spatial anisotropy of 
the orbitals, it is found that the state with the d xz ferro-orbital order and AFM spin 
order will have the lowest energy. The same argument can be applied to NNs along 
the y direction. Along the diagonal direction between next nearest neighbors (NNN) , 
the inter-orbital hopping no longer vanishes and there is no dominant hopping 
amplitude. Therefore we may make the simplification that the exchange constant is 
orbitally independent and equals J2 no matter which orbitals are occupied. 
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Fig. 3. Orbital order and the resulting magnetic order with the anisotropic exchange constants, 
(a) At high temperature, d xz and d yz orbitals are equally populated, (b) The d xz orbital order 
and the Q = (tt,0) AFM order, (c) The d yz orbital order and the Q = (0,7r) AFM order. J la , 
Jit, and J2 are the exchange constants between NNs and NNNs. 



As a result, we can write down the following spin-orbital model to desc ribe the 
structural-magnetic phase transitions in iron-based superconductors) ^^ 



Jla 



ff + - 

1 2 



i-\-x 



Jib [T?-- 



E 

i 

E 



Jib t. 



i+y 



Jla T, 



Si • Si- 



Si ■ Si- 



JiSi ■ Sj 



(9) 



We have ignored the weak FM exchange from Eq. ([8]), and the orbital pseudospin 
is reduced to an Ising variable. At high temperature, the orbital is disordered, 
with each sites being equally populated by d xz and d yz , as illustrated in Fig. [3ja). 
Consequently, the magnetic exchange is weak and there is no magnetic order. At low 
temperature, we have a ferro-orbital order of either d xz or d yz [Fig. [3fb) and (c)]. 
The effective spin model becomes an anisotropic Ji a -Jib-J2 Heisenberg model. If 
J2 > Jib, we have a stripe- like AFM order with the ordering wave vector Q = (ir, 0) 
or (0,tt). 

We note that our analysis is based on strong anisotropy in the hopping am- 
plitudes of the d xz and d yz orbitals. This assumption is confirmed by DFT 



calculations^ i n which a similar scenario is also proposed. Experimentally inelas- 
tic neutron scattering resulti^ give further support that the spin physics in iron 
pnictides is governed by a strongly anisotropic Heisenberg model. Now the question 
is whether orbital order can exist in the absence of the long-range magnetic order. 
First, we note that the orbital order is an Ising order whereas the spin order is 
an 0(3) order. In a true two-dimensional system, orbital order can occur at finite 
temperature whereas spin order only sets in at T — 0. Turning on couplings along 
the third dimension, we can drive the magnetic transition close to the orbital tran- 
sition. Nevertheless, there exists the nematic phase with nonzero orbital order but 
no magnetic order. This argument is also used by some theories on the basis of 
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Fig. 4. Illustration of the Fermi surface distortion from orbital order, (a) Fermi surfaces without 
orbital order, (b) Fermi surfaces with d xz ferro-orbital order, (c) Fermi surfaces with d yz ferro- 
orbital order. In (b) and (c) the original Fermi surfaces are shown by the dashed lines. The arrows 
represent the nesting instability between the hole and electron pockets. The color represents the 
dominant orbital at that location, with 'red' corresponding to d X z and 'blue' to d yz . 



spin fluctuations!^^ Numerically, such a sequence of phase transitions is found in 
both exact diagonalizatiorf^ and Monte-CarlcJ^ studies of the spin-orbital model. 
Second, from our discussion, orbital order can be driven by other interactions, like 
electron-phonon and quadrupole-quadrupole interactions, in addition to the spin- 
orbital couplings. Therefore orbital order transition can have a distinct energy scale 
from the magnetic transition, and happen at a much hig her temperature. This is 
indeed the case in many other transition metal oxidesl^l 



2.3. Orbital order from weak coupling 

In the previous subsection, we have described how orbital order arises in the strong- 
coupling limit. However, the normal state of the iron-based superconductors remains 
metallic. As a consequence, the opposite weak-coupling limit is worthy of being 
studied as well. 

Let us examine the Fermi surface of the iron-based superconductors. In the 
unfolded Brillouin zone that contains one Fe atom per unit cell, there is a circular 
hole pocket at the zone center T = (0,0), whereas two elliptical electron pockets 
appear at X = (n, 0) and Y = (0, 7r), respectively [see Fig.@|a)]. Here for simplicity 
we have ignored the second hole pocket at T and a third possible hole pocket at the 
zone corner M = (ir, ir). In terms of a two-orbital model, the eigenstates on different 
locations of the Fermi surfaces have different weights of d xz and d yz orbitals. One 
prominent feature here is that under a tt/2 rotation, the shape of the Fermi surface 
is invariant, but the dominant orbital weight changes between d xz and d yz . This 
condition is strictly enforced because the two orbitals exchange with each other 
under the C4 rotation. 

It is then obvio us that an electron nematic order can be captured by a ferro- 
orbital order ESlQ] a s shown in Fig. HJb), if d xz orbital has a larger occupation 
number, the hole pocket will be elongated along the x direction and shortened 
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along the y direction. Such a Fermi surface distortion is a hallmark of the nematic 
order. Meanwhile, the electron pocket at X shrinks while the one at Y enlarges. 
Similarly for the d yz orbital order, the opposite takes place [Fig. QJc)]. So in this 
multiorbital system, orbital order naturally leads to an electronic nematicity. 

This ferro-orbital order has an important consequence in picking up the preferred 
wavevector of the stripe-like AFM order. In this itinerant model, the magnetic order 
is induced by the nesting between the hole and electron pockets. As we can see in 
Fig. EJa) , for the isotropic Fermi surface, the nesting can happen between either of 
the two electron pockets. Furthermore, due to the mismatch between the hole and 
electron pockets, the nesting condition is far from perfect and the nesting wavevector 
is usually incommensurate. In the presence of the ferro-orbital order, the distorted 
Fermi surface enhances the nesting condition along one of the two directions with a 
commensurate nesting wavevector. In particular, d xz ferro-orbital order favors the 
(7r, 0) magnetic order while the (0, n) order arises in the d yz orbitally ordered state. 
This argument also gives a natural explanation why the structural phase transition 
always precedes the magnetic transition. 

Very recently, Avci et. al. observed a novel first order transition belo w T /y into 
a new magnetic order with C4 symmetry restored in Bao.76Nao.24Fe2As2!^ While 
the exact configuration of this new magnetic state remains to be determined, it has 
been shown that a spin-nematic model could lead to several candidate states within 
a limited range of parameter space. We note that such 're-entry' of C4 symmetry 
could still be consistent with an orbital model. Let us denote M\ as the stripe- 
like AFM order with orbital order and M2 as this new magnetic state with C4 
symmetry Since M\ and Mi break different symmetries, a first order transition 
between them could naturally occur if they are competing orders. This is a general 
principle regardless the nature of the model. As a result, if a spin-nematic model 
could produce such a first order transition, so could an orbital model in principle. 
In fact, this again reflects that the mean- field phase diagram is insentitive to the 
details of the model as long as the symmetries of the phases are captured correctly. 
Nevertheless, it is still intriguing to further explore the physical properties of this 
new magnetic state both experimentally and theoretically. 

It is then clear that in order to resovle the microscopic mechanism driving the 
phase transitions in iron-based superconductors, the investigation on quantum fluc- 
tuations beyond the mean-field states is necessary. As will be shown in the next 
section, the ferro-orbital order in a multiorbit al sys tem is in fact identical to the 
Pomeranchuk instability in the d-wave channel! ^^ The elementary excitations in 
the proximity of the orbital ordering quantum critical point are unusual, leading to 
a non-Fermi liquid behavior which cannot be caused by the spin nematicity. 
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3. Excitations near the orbital ordering quantum critical point 

3.1. Non- Fermi liquid behavior near quantum critical point in 
itinerant systems 

In this subsection, we briefly review the principle behind the non-Fermi liquid be- 
havior emerging from the quantum criticality, which has been a hot research topic 
over the last decades. Generally speaking, a phase transition can be described by a 
problem of minimization of the free energy F = U — TS, where U is the internal 
energy of the system, T the temperature, and S is the entropy of the system. In clas- 
sical systems, the thermal fluctuations are dominant and hence, the phase transition 
is a compromise between states with lower U and higher S. At high temperature, 
F is minimized by the state with higher entropy which has the most symmetries 
the system can have. On the other hand, at low temperature a state with lower 
internal energy is favored. In an interacting system, such a state could have less 
symmetries. As a result, a phase transition occurs when the system spontaneously 
loses some symmetries at a critical temperature T c as it is cooled down from high to 
low temperature. An indicator for the phase transition is an order parameter which 
develops a non-zero expectation value only in the low temperature ordered state. 

In real systems, we can have a tuning parameter, namely <?, that we can use 
to change T c . Typical examples for the tuning parameter are ext ernal pressure, 
magnetic field, and chemical dopings. A quantum phase transitiorP31 is then referred 
to a transition occurring at a quantum critical point (QCP) defined at T c — and 
g = g c . A surprising feature of this subject is that although the quantum phase 
transition is strictly at zero temperature, the associated quantum fluctuations near 
QCP can fundamentally change the physical properties even at finite temperature. 
This is why the field of quantum criticality has been so vibrant recently. Along with 
its theoretical novelty, it is also very relevant to the strongly correlated materials. 
A recent review on new progress in several QCPs of insulating phases can be found 
in Ref. HO 

While the quantum fluctuations associated with the QCP can be described by 
the dynamics of the order parameters, the situation is further complicated in an 
itinerant system. In an insulating phase, by definition, all the quasiparticle spectrum 
is gapped, and as a consequence, the low energy physics contains only the collective 
excitations from the order parameter dynamics. If the QCP, however, resides in an 
itinerant system, the existence of the gapless quasiparticles can affect the dynamics 
of the order parameter and in feedback, the quasiparticle properties can also be 
altered fundamentally. The pioneering wor k in demonstrating this idea was done by 
HertzS^l and later improved upon by Millis^l 

Consider a Fermi liquid near a ferromagnetic QCP. Since the ferromagnetic order 
is at zero momentum (q = 0), we analyze the low energy excitations near q = 0. 
The particle-hole continuum from quasiparticles is bounded below Max\e(k + q) — 
e(fc)] ~ \q\, where e(fc) is the quasiparticle energy at momentum k. Given that 
the dispersion of the magnon in a ferromagnet goes like q 2 , as the ferromagnetic 
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QCP is approached from the Fermi liquid side, a branch of spin collective modes 
emerges with a dispersion of 5(g) + aq 2 , where 5(g) can be viewed as the mass of the 
magnon that vanishes at the QCP (5(g —> g c ) —> 0). Clearly, right at the QCP the 
whole branch of the spin excitation modes lies completely inside the particle-hole 
continuum and hence, these collective modes can no longer have infinite life time 
and will eventually decay into uncorrelated particle-hole pairs. As shown by Hertz 
and Millis, this overdamping effect brings in a correction to the dynamical exponent 
z and can be captured even at the random-phase approximation (RPA) level. 

Now we are ready to study the fate of the quasiparticles. Since these overdamped 
spin modes near zero momentum still have non-zero spectral weights, quasiparticles 
can be scattered off them. As a result, the life time of the quasiparticles, which is 

inversely proportional to the imaginary part of the quasiparticle self-energy ImS(w), 

J47I 

is no longer simply determined by the phase space argument as for a Fermi liquid . 
Rather, it is mostly determined by the dimension of the system as well as z of the 
overdamped collective modes, ImE(w) ~ w^ 2 , and in the current two-dimensional 
ferromagnetic QCP case, ImE(u) ~ w 2 / 3 . This is the milestone result signaling 
the violation of the Fermi liquid behavior (ImS(w) ~ lo 2 ) induced by the quantum 
criticality. 

Although the example given above is for a ferromagnetic QCP, the same the- 
ory can be applied to all the zero-momentum ordered states transitioned into 
from a Fermi liquid since they fall into the same effective field theory in Landau- 
Ginzburg formalism. Another exampl e is the two-dimensional nematic quantum 
fluid fi rst proposed by Oganesyan et al.\^ which motivated a number of theoretical 
efforts 48 49 5 °EIEI As discussed in Sec. [H the nematic order is a translationally 
invariant phase breaking the rotational symmetry in which the Fermi surface is spon- 
taneously distorted along one particular direction in the two-dimensional plane. In 
fact, the instability toward the change of the Fermi surface was studied a long time 
ago in 1958. PomeranchulP2l first derived the criteria for the Fermi surface instabil- 
ities induced by the interactions at higher angular momentum channels described 
by the Landau parameters F n . The nematic order is one of these Pomeranchuk in- 
stabilities which is in the d-wave interaction channel F2 , and it has received a lot of 
attentio n in the la st decade due to the discoveries of anisotropic phases in several 
systems!21^2154155| What was unexpected back at the time of Pomeranchuk is that 
the nematic quantum fluid could exhibit a non-Fermi liquid behavior in the same 
way as the itinerant ferromagnet, except that the interaction in this case would be 
from the charge instead of the spin. The reasoning is as follows. Because the ne- 
matic order still preserves the translational invariance, it must be an order at zero 
momentum. Moreover, since it breaks a continuous symmetry of rotation, it will 
develop Goldstone modes with a dispersion of q 2 . Finally, since the Fermi surface is 
only distorted but not killed, there are still itinerant electrons and consequently the 
particle-hole continuum is still present. As a result, all the physics discussed in the 
ferromagnetic QCP can be applied to the nematic QCP despite some differences in 
the details. Indeed, the overdamped collective modes with z — 3 as well as the non- 
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Fermi liquid behavior were found by Oganesyan et alW^ using RPA theory, and later 
Lawler et al. confirmed these findings using higher dimensional bosonization!^ We 
emphasize that the consistency between the results from RPA and higher dimen- 
sional bosonization is remarkable since the former is perturbative while the latter 
is non-perturbative. This suggests that RPA theory is asymptotically exact at long 
wavelength, and an analysis based on RPA theory for the quantum criticality near 
zero momentum is a reasonable starting point despite its perturbative nature. 

There is a concern that since the materials of interest are crystalline and do not 
have continuous rotational symmetry, the Goldstone modes are not massless, which 
would be at odds with the occurrence of the non-Fermi liquid. The main effect of 
this is that at zero temperature, in a lattice model the non-Fermi liquid behavior 
goes away immediately in the nematic ordered phase, while in a continuous model 
it survives throughout the ordered phase. However, as long as the transition is of 
a second order, the critical fluctuations are still overdamped and massless at the 
QCP in a lattice model, and the non-Fermi liquid behavior still prevails. At finite 
temperature, if the thermal fluctuations are large enough, i.e., fcsT c > m where m 
is the mass of the Goldstone mode, the non- Fermi liquid behavior still emerges. The 
non-Fermi liquid behavior near t he nem atic QCP at zero and finite temperature in 
a lattice model has been 

studiedPSMl 
3.2. Analysis on a two-orbital model 

As discussed in the previous sub-section, the finger print of the non-Fermi liquid 
behavior is the collective mode arising from critical fluctuations near the QCP with 
a dynamical exponent of z = 3. This is commonly the case for a zero- momentum 
ordered state. In this section, we will review the study on a two-orbital model which 
is the minimal model for the iron pnictides and bilayer Sr3Ru207. The key point 
here is that the orbital order in a two-orbital model is identical to the nematic 
order in a single-band lattice model. Hence they share the same physics, including 
the non-Fermi liquid behavior near the QCP. 

To see this connection, we start from the Fermiology of the two-orbital model 
discussed in Sec. 2. Within the framework of the Fermi liquid theory, the interactions 
that are relevant should be the ones between states on the Fermi surface. Without 
a loss of generality, th e eigenstate wavefunction on the a Fermi surface can be 
generally written as!^SI 

\ip a ,a(p)) = (coscf>p\d xz (ft)) +sixi(f>p\dy Z (p))) ®x<j (10) 

where <j)p is the hybridization angle between d xz (p) and d yz (p) at p and \ a is spin 
eigenstate. In real space, the C4 symmetry is enforced by the interchange of d xz and 
d yz orbitals up to a sign under a -| rotation, R^ . This places a relation between 
the eigenstate wavefunctions at p and p 1 — R^.p of = (j>p + ^. Following the 
standard procedure, we evaluate the interaction V(pi,p>2) between states at p\ and 
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p 2 on the Fermi surface as 



(i) a ,a(Pl),^ a ,*(P2)\V(p 1 ,p 2 )\i) a ,a(P2),i> a ,a(Pl)}- 



(11) 



Eq. (|llj) can be evaluated by the Hartree-Fock approximation which gives 



Clearly, the Hartree term is not affected by the multiorbital nature, but the effect of 
the orbital hybridization comes in the Fock exchange part, leading to the effective 
Landau interaction function, 



Now it becomes apparent that because of the form factor cos2(^j —<t>p 2 ), the orbital 
hybridization induces a d-wave Landau interaction regardless the original form of 
V(p\,p 2 ). In other words, the tendency toward the nematic order, or equivalently 
an orbital order, is naturally enhanced in such a two-orbital system. 

Applying the same procedure outlined above, one can readily show that the 
two-orbital model can be mapped to a two-band model with a large F 2 interaction 
on its Fermi surfaces. Such a mapping has to be true because orbital degrees of 
freedom are only well-defined on a lattice. As a result, it will eventually evolve 
into a band model when approaching the continuum limit. Since there is no phase 
transition induced simply by taking the continuum limit, the orbital model has the 
same physics as its corresponding band model. This indicates that the orbital order 
in a weak-coupling itinerant system is identical to the nematic order from the d- 
wave Pomeranchuk instability. This is different from the strong coupling limit in 
which the orbital order is an Ising-like order as discussed in Sec. 12.21 This should 
not be surprising since quite generally the hybridizations between different orbitals 
are in the single particle part of the Hamiltonian, not in the interaction part. As a 
result, the orbital index is not a good 'quantum number' in the weak coupling limit 
but remains approximately a good one in the strong coupling limit. 

There have been several studies on the two-orbital model beyond the mean-field 
theory. A renormalization group analysis showed that the orbital order instability 
can be largely enhanced by spin fluctuations, which has bee n expl oited to ex plain 
the magnetic field-induced nematicity observed in S^Po^O?! ^^ Lo et aZl^have 
analyzed the two-orbital model using higher dimensional bosonization. Near the 
orbital ordering QCP, it is found that a branch of overdamped collective modes 
with z — 3 does emerge. This provides a strong foundation for the occurrence of 
the non-Fermi liquid behavior in the iron pnictides which will be reviewed in the 
next section. 




V{pi 7 p 2 ), 

~V(p 1 ,p 2 )\(lp a ,a{Pl)\lp a ,ai.P2))\ 2 - 



(12) 



Fifufc) = V(q = 0) - i[l + cos2(0 p - - fa)} x V(p u p 2 ), 
f a (Pi,P2) = -i[l + cos2(^ -fa)} x V(p t ,n). 



(13) 
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3.3. RPA study on a five-orbital model 

Iron pnictides are known for their complicated multiorbital structure. It has been 
suggested that in the unfolded Brillouin zone, a model containing five orbitals is a 
much better choice to capture the details in the electronic structurePZl Although 
the orbital order involves mostly the d xz and d yz orbitals, the couplings to other 
orbitals can substantially modify the hybridization between d xz and d yz . Therefore 
a theory based on a five-orbital model capturing enough detail is certainly necessary 
to demonstrate the non-Fermi liquid behavior. Fortunately, since the RPA theory 
is a reliable starting point as explained in Sec. 3.1, one can directly calculate the 
single particle self energy^l a t the RPA level with a multiorbital Hubbard modelPSl 
The computational scheme for the self energy Y, a t(k,ui) can be depicted by the 
Feynman diagram plotted in Fig. El and we refer the readers to Ref. 60 for de- 
tails. It is showrP^l that ImS n i,(^if,w) changes from a Fermi liquid (~ uj 2 ) to a 
non-Fermi liquid (~ uj a ,a < 1) behavior as the orbital ordering quantum critical 
point is approached within a realistic range of interaction parameters for the multi- 
orbital Hubbard model. This further confirms that the important physics from the 
two-orbital model survives in a complicated five-orbital model. In Sec. 21 we will 
discuss that such a non-Fermi liquid behavior gives a natural explanation to the 
novel enhanc eme nt of the zero bias conductance observed in a recent point contact 
measurement 



It should be noted that in this five-orbital Hubbard model, contributions to the 
single particle self energy from the spin interaction are actually not small. As a 
result, it is an open debate as to whether the orbital ordering is caused from charge 
or spin interactions, and with the current model, there seems to be no hope in 
resolving this issue. Nevertheless, if one includes the nearest-neighbor interaction, 
the contributions from charge interaction can be enhanced. This is demonstrated^ 
within the context of Sr3Ru207 which is another intriguing material exhibiting 
the anisotropic phase under magnetic field. Moreover, Kontani et al. showed that 
including the electron-phonon coupling as well as the vertex corrections beyond the 
RPA level can boost up the orbital orde r instabili ties in the charge channel which 
may lead to an even richer phase diagramJ 63 | 64 | 65 [ gj nce these new additions to the 
multiorbital Hubbard model are all realistic and not negligible in general, the orbital 
ordering from charge interaction should be more important than it now appears, at 
least for some families of the iron pnictides. 

The non-Fermi liquid behavior could be a strong indicator for the orbital fluctu- 
ations because fluctuations from spin nematicity do not seem to produce the same 
result. The reason is given as follows. As described in Sec. [2j in the spin nematic 
scenario the nematic phase breaks the symmetry arising from the degeneracy 
of two possible ordering wavevectors (ir,0) and (0, ir). Following this theory, the 
critical fluctuations near the nematic QCP would mostly arise from the spin flue- 
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Fig. 5. Feynman diagrams for the single particle self energy at the RPA level. The dashed line 
refers to the interactions between different orbitals, and the solid line is the single particle propa- 
gator G aa . All the symbols refer to an index of orbital and spin. The Greek symbols are internal 
indices which are summed over all orbitals and spins. 

tuations around finite momenta (71", 0) and (0, 7r), which, according to Hertz-Millis 
theory, usually have a much weaker capability of inducing non- Fermi liquid behav- 
ior. To date, we are not aware of any theoretical study proposing a non- Fermi liquid 
behavior arising from spin nematicity. 

4. Experimental signatures for orbital fluctuations 
4.1. Point contact spectroscopy 
4.1.1. Introduction 

A point contact is simply a contact between two metals whose characteristic size 
d is much less than the electron elastic and inelastic mean free paths: l e i, U n ^> d. 
Point contact spectroscopy (PCS) studies the non-linearities of the current- voltage 
(I-V) characteristics of these metallic constrictions. The bias voltage applied to the 
point contact junction determines the energy scale of the scattering process. If both 
metals are ohmic, the I-V curve is linear and the differential conductance dl/dV 
is constant. However, when scattering processes are present, they show up in t he 
PCS spectrum. PCS was initially used to study electron-phonon interactions^^ in 
metals. If one of the metals is replaced by a superconductor, making a N-S junction, 
Andreev reflection^! is observed. In recent years, PCS has also been used to probe 
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heavy fermion compounds where i t is sensitive to the onset of the Kondo lattice 
that appears as a F ano line shape^l and to the hybridization gap in the heavy 
fermion URu2Si2.^ PCS has also been used extensively to study the iron based 
superconductors and provides evidence for orbital fluctuations in the normal sate 
of underdoped Ba(Fei_ x Co x )2As 2 , parent SrFe2As 2 and Fei +y TeP^ 



4.1.2. Conductance Spectra 

Arham et alW^ reported conductance spectra on Ba(Fei_ x Co x ) 2 As2, spanning the 
entire phase diagram. Part of their data is reproduced in Fig. [51 For the undoped 
parent compound [Fig.[6ja)], at the lowest temperature (blue curve), they saw a dip 
at zero bias and two asymmetric conductance peaks at ~ 65 mV. This double peak 
feature was superimposed on a parabolic background. As the temperature was in- 
creased, the dip at zero bias filled up, the conductance peaks moved inward, and the 
bias voltage range of the conductance enhancement decreased. No dramatic change 
occurred as the structural transition temperature, Ts, was crossed (red curve). The 
enhancement eventually disappeared at 177 K, more than 40 K above Tg. 

For underdoped Ba(Fei_ x Co x )2As2 [Fig. E[b)], where superconductivity coex- 
ists with long-range magnetic order, they observed Andreev reflection at low voltage 
biases in the superconducting state. However, just like the parent compound, two 
conductance peaks occurred at ~ 65 mV. Above the onset temperature of the super- 
conducting transition T c , Andreev reflection completely disappeared and the high 
bias conductance evolved just as it did for BaFe2As2- The right inset in Fig. [S^b) 
shows a zoom in of the Andreev reflection features while the left inset plots the 
conductance spectra above T c on a log plot. 

For overdoped Ba(Fei_ x Co x )2As2 [Fig. EJc)], Andreev reflection was observed 
in the superconducting state, but unlike the underdoped compounds, no higher 
bias conductance peaks were detected. Above T c , a parabolic background remained 
which flattened with further increase in temperature. 

Conductance spectra were also presented on hole underdoped Bai_ x K x Fe2As2 
[Fig. |UJd)]. The sample has a coexistence of magnetism and superconductivity (Tjv 
= Ts ~ 90 K, T c ~ 20 K). Below T c clear Andreev reflection was observed. Above 
T c , Andreev reflection disappeared, leaving a downward facing background that did 
not change with any further increase in temperature. This is remarkably different 
from the situation in electron underdoped Ba(Fei_ x Co x )2As2. 

As comparisons, CaFe2As2 and SrFe2As2 were also probed. The data for 
SrFe2As2 [Fig. 2(a) in Ref. IBTj were very similar to those of BaFe2As2, with conduc- 
tance enhancement around zero bias that sets in before the magnetostructural tran- 
sition. However, for CaFe2As2 [Fig. 2(b) in Ref. 161). the conductance enhancement 
disappeared around 100-110 K, which is much lower than the magnetostructural 
transition temperature, 170 K. 

Similar features were also observed in the parent compound for the iron chalco- 
genide superconductor. Fei.i3Te showed a conductance enhancement that survived 
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Fig. 6. From Ref. 61, copyright 2012 by the American Physical Society, (a) Conductance spectra 
for BaFe2As2. Conductance enhancement with peaks at ~ 65 mV superimposed on a parabolic 
background was observed at low temperatures. The peaks moved in as the temperature was in- 
creased and the enhancement survived well above Ts (red curve), (b) Ba(Feo.gsCoo.o5)2As2 dis- 
plays a coexistence of magnetism and superconductivity. At low temperatures, clear Andreev 
peaks were observed [right inset (b); the arrows are pointing out the Andreev peaks]. A conduc- 
tance enhancement with peaks at ~ 65 mV coexisted with the Andreev spectra and evolved with 
temperature as it did for BaFe2As2. This enhancement increased logarithmically near zero bias 
[left inset (b)]. The overdoped compound Ba(Feo.875Coo.i25)2As2 showed Andreev spectra below 
T c . It did not have conductance peaks at higher bias values like the Co underdoped compounds, 
(d) The hole underdoped Bao.sKo.2Fe2As2 has a coexistence of superconductivity and magnetism. 
It showed Andreev spectra below T c and no higher bias conductance enhancement. This was in 
contrast to the data obtained from electron underdoped Ba(Fei_ x Co x )2As2 [F'ig. [6fb)] . 



above the magnetic and structural transition temperatures [Fig. 2(c) in Ref. W\\ . 
The conductance enhancement was observed till 75 K (T/v = Ts ~ 59 K). 

To summarize the work of Arham et ai.P^ BaFe2As2, SrFe2As2, underdoped 
Ba(Fei_ x Co x )2As2 and Fei+ y Te exhibited a conductance enhancement that sets in 
above Ts, CaFe2As2 only showed the enhancement below Ts while Bao.sKo.2Fe2As2 
did not show any conductance enhancement. Overdoped Ba(Fei_ x Co x ) 2 As 2 does 
not have a Ts and only showed Andreev spectra below T c . The high bias background 
for all compounds except for Bao.sKo.2Fe2As2 was an upward facing parabola. 
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4.1.3. Discussion 

Arham et aZ.^J noticed a correlation between the presence of conductance enhance- 
ment around zero bias and in-plane resistivity anisotropy in the compounds. For 
detwinned under doped AEFe 2 As2 it has been shown that below Ts a resistivity 
anisotropy exists. | 10|H|16 170| Above Tg there is notable anisotropy for AE = Ba, 
negligible anisotropy for AE = Sr, and no anisotropy for AE — Ca. Detwinned 
Fei_|_ y Te also shows a resistivity anisotropy above the structural transition!-^ The 
anisotropy above Tg is sensitive to the uniaxial force required to detwin the sam- 
ples. Detwinned underdoped Bai_ x K x Fe2As2 does not show any anisotropy at all, 
cither below or above TgM^ 1 

The presence or absence of the in-plane resistivity anisotropy matches whether 
or not a conductance enhancement is detected. The correlation of the conductance 
enhancement with the resistivity anisotropy indicates they are likely caused by 
the same underlying physics. Arham et al.-- constructed a revised phase diagram 
for Ba(Fei_ x Co x )2As2 marking a new line on the underdoped side showing the 
temperature below which the co ndu ctance enhancement was observed [Fig. 0a)]. 

Theoretical work by Lee et alW^ (reviewed in Sec. 3) showed that orbital fluctu- 
ations above Tg were expected to provide extra contributions to the single-particle 
density of states (DOS) at zero energyp22 The DOS followed a log dependence as 
the energy was increased. Arham et al. compared their data with this prediction 
and found that the conductance enhancement for BaFe2As2 above Tg followed a 
log dependence from ~ 40 mV to ~ 90 mV [Fig. Eft))]. They obtained similar fits 
above Tg for SrFe2As2 and Fei.i3Te. Furthermore, the absence of similar effects in 
the data on Bao.sKo.2Fe2As2 was consistent with the prediction that crystals that 
did not show the resistivity anisotropy would also not exhibit the excess conduc- 
tance due to orbital fluctuations. Their data therefore strongly indicated that the 
enhancement in conductance observed was a consequence of orbital fluctuations. 

It should be kept in mind that the conductance (dl/dV) measured by point- 
contact spectroscopy does not directly correspond to the density of states. PCS 
data is a convolution of the Fermi velocity and the energy-dependent density of 
states along with any scattering processes that might be present. For normal metals, 
the Fermi ve locity and the density of states are inversely related and cancel each 
other outP^ There is a lack of theoretical models for interpreting PCS data on 
correlated metals, where the DOS are energy dependent and do not cancel out with 
the Fermi velocity when dl/dV is measured. A theory considering both the energy 
dependence of the electronic DOS and scattering processes would be extremely 
helpful in obtaining a better understanding of the experimental data. 
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Fig. 7. From Ref. |61l copyright 2012 by the American Physical Society, (a) Phase diagram for 
Ba(Fei_ x Co x )2As2 marking a new line on the underdoped side showing the temperature be- 
low which the conductance enhancement was observed by PCS. (b) Conductance above Tg for 
BaFe2As2 followed a log dependence from ~ 40 mV to ~ 90 mV. 



4.2. Magnetic torque measurement on BaFe2(Asi_ x P x )2 
4.2.1. Introduction 

Torque magnetometry is an excellent method to detect the presence of C4 symmetry 
breaking novel electron matter in a tetragonal crystal lattice. The quantity to be 
measured is the torque defined as f = uqVM x H, where no is the permeability of 
vacuum, V the crystal volume, and M is the magnetization induced by magn etic 
field H. When H is rotated within the tetragonal a-b plane, it can be showrff^ 
that the torque has a component of T20 = \noH 2 V[(xaa — Xbb) sin 2(f) — 2\ a b cos 20], 
where <p is the azimuthal angle measured from the a axis and Xij i s determined 
by Mi — J2jXijHj- When the system has tetragonal symmetry, \ aa — \bb and 
Xab = 0, making r 2 = 0. A non-zero r 2( p, meaning either Xaa ^ Xbb or Xab 0, 
signals tetragonal (C4) symmetry breaking. 



4.2.2. Experimental data and discussion 

Kasahara et ali^ carried out torque magnetometry measurements on 
BaFe2(Asi_ 2; P a ;)2 spanning the entire phase diagram. They observed that ob- 
tained a non-zero value at an onset temperature T* higher than Tg, indicating 
a nematic phase transition occurring above the structural phase transition. They 
further showed that the onset temperature, T*, of this new nematic phase was sup- 
pressed with increasing P doping but persisted well into the overdoped side. Above 
T* the torque signal is isotropic. 

The fact that T20 becomes non-zero before the structural transition temperat ure 
is a clear indication of the presence of an electron nematic phase. Kasahara et oZP^ 
noticed that below T*, the function form of the was: = A 2( p cos 2ip. This 
meant that Xaa = Xbb, Xab 7^ 0, and the nematicity was along the [110] (Fe-Fe 
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bond) direction. 

Kasahara et al. also analyzed x-ray data on these crystals using a two-peak 
fitting procedure (Fig. 3 in Ref. [T7|). They observed a very small but finite lattice 
distortion in the temperature range T$ < T < T*. They argued that this was 
due to a linear coupling between the order parameter governing the orthorhombic 
lattice distortion and the one representing electronic nematicity. They postulate 
that the true thermodynamic transition thus occurs at T* and the one at T$ is a 
meta-nematic transition. 

The origin of the electronic nematicity occurring at T* is still under debate. The 
magnetic field applied during the torque magnetometry measurement breaks the C4 
symmetry itself since it points along one direction on the two-dimensional plane. 
Thus there is an uncertainty whether the transition occurring at T* is a true second 
order phase transition or not. Nevertheless, there is no doubt that certain fluctu- 
ations associated with the nematicity must be present above the structural phase 
transition. While spin fluctuation is a plausible candidate for causing nematicity in 
the underdoped samples, the observation of nematicity in the region far from the 
antiferromagnetic phase casts a strong doubt on this. Keeping this in mind, orbital 
fluctuations might be a better candidate for causing the nematic phase. In particu- 
lar, it is intriguing to see that T* for the parent compound BaFe2As2 is remarkably 
consistent with the onset temperature for the enhancement of the zero-bias conduc- 
tance discussed in Sec. 4.1. As a result, in the orbital-based scenario, T* would be 
the temperature for a crossover at which the orbital fluctuations become important, 
instead of the transition temperature for a second order phase transition. 

4.3. Inelastic neutron scattering measurement on 
Fe 1+y _ x (Ni/Cu) x Te . 5 Se . 5 

4.3.1. Introduction 

Inelastic neutron scattering measurement on unconventional superconductors pro- 
vides information on their bosonic (phonons and magnons) excitation spectra. After 
intensive international efforts in applying this technique to cuprates, heavy fermions, 
and iron based superconductors, the behavior of the spin excitations for magnetic 
fluctuation driven superconductors has been well established. In the normal state, 
the low energy spin fluctuations usually have a peak feature near a certain wavevec- 
tor which is typically the ordering wavevector in the parent compound, e.g., (tt,tt) 
for cuprates and (77, 0) for the iron based superconductors. The appearance of this 
peak is an indicator of the strong spin fluctuations inherited from the magnetic or- 
der in the parent compounds. On cooling below T c , the low energy spin fluctuations 
at the same wavevector develop a gap due to the superconductivity, and the spec- 
tral weight gets shifted to a pronounced resonance peak. However, qualitatively, the 
magnetic spectrum above and below T c tends to be the same. Moreover, if there is 
no other phase transition in the normal state, the magnetic spectrum is unchanged 
with rising temperature, and the only prominent effect is the broadening of the 
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spectral weights due to the thermal fluctuations. A nice review can be found in 
Ref.ffl 

4.3.2. Experimental Data 

Xu et a?P31 probed the low energy spin fluctuations in Fei+ y _ x (Ni/Cu) x Teo.5Seo.5 
using inelastic neutron scattering. The Ni/Cu substitution of Fe was done to perturb 
the system and reduce T c . The unperturbed system has a T c of 14 K, 2% Ni brings 
it to 12 K, 4% Ni to 8 K, while 10% Cu suppresses superconductivity entirely. 
The magnetic neutron scattering intensity was probed near the M point, Qaf = 
(0.5,0.5,0), along with its evolution along the wave vector q = [1,-1,0]. It was 
found that the temperature dependence of the magnetic excitation spectrum was 
qualitatively different for superconducting and non-superconducting samples. In the 
superconducting sample (4% Ni doping), the low energy spin excitation below T c 
showed a commensurate U-shaped dispersion, with the bottom of the U close to the 
spin resonance energy mode. As expected, the only difference on crossing T c was the 
change in the spectral weight of the resonance peak. However, at 100 K, the spectra 
looked remarkably different, showing two columns centered around incommensurate 
wavevectors Q = (0.5 ± 5, 0.5 =F 6,0) at low energies. Xu et al. described it as 
looking like the 'legs of a pair of trousers'. For the non-superconducting sample 
(10% Cu doping), the trousers shaped spectra were obtained at both low and high 
temperatures. In other words, the prerequisite condition for superconductivity in 
this family is the occurrence of an incommensurate-to-commensurate transformation 
in the magnetic excitation spectrum in the normal state. For the superconducting 
samples, this crossover appeared to occur close to 3T C . 

4.3.3. Discussion 

The question confronting Xu et al. was how the low energy magnetic spectra could 
evolve from two incommensurate vertical columns at T > T c to a commensurate 
U-shaped dispersion at low temperatures. If one considers only the spin interac- 
tions, such a transformation seems to be impossible without a magnetic transition. 
However, it can be naturally understood if there exists another energy scale which is 
distinct from the spin sector but could affect the spin excitation. As shown in Sec. 2, 
in an itinerant system orbital order can change the Fermi surface. Therefore the spin 
excitation spectrum, which is basically the flip of the electron spin as it is scattered 
between states on different parts of the Fermi surface, can be changed qualitatively 
by th e or bital order and related fluctuations. This has been demonstrated by Lee 
et al.\^ and the incommensurate-to-commensurate transformation is reasonably 
reproduced from a five-orbital model including Gaussian fluctuations to account 
for the low energy fluctuating orbital correlations. This orbital correlation-induced 
transformation can be viewed as an analogy of the stripe-induced neutron anomaly 
in cuprates, in which the stripe order or correlation renders the neutron peak from 
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(tt, 7r) to incommensurate wavevectors. 
5. Concluding Remarks 

In this paper, we have reviewed some topics of the orbital physics originating from 
degenerate d xz and d yz orbitals. Although the current focus is on understanding 
some non-trivial experimental results in the iron pnictides, these theories are gen- 
erally applicable to all systems having degenerate d xz and d yz orbitals on the Fermi 
surfaces, e.g., Sr3Ru20y. While we fully recognize the importance of spin physics in 
the iron pnictides, our purpose is to point out that there exist some features and ex- 
perimental data which seem to be understood more naturally by the orbital physics. 
In particular, since iron pnictides are generally less correlated compared to cuprates, 
no energy scale is really dominant over any other. This strongly suggests that many 
physical properties are probably not universal and as a consequence, it is highly 
possible that some families of the iron pnictides have strong orbital physics charac- 
ter while others do not, depending on their detailed electronic structures. To achieve 
a better understanding, more spectroscopic measurements are certainly necessary. 
For example, a recent Raman measurement has shown a superconductivity-induced 
peak in the Bi g channel in Bao.6Ko.4Fe2As2, which was originally interpreted as the 
existence of a competing d-wave superconducting component!^ One of us, however, 
showed that the orbital fluctuations in the superconducting state could lead to such 
a peak as wellP^ Nevertheless, it appears to be surprising that superconductivity 
is ubiquitous among the iron pnictides despite this non-universality. If it turns out 
that superconductivity does arise from a normal state with strong orbital fluctua- 
tions in some families of the iron based superconductors, a unified picture for the 
pairing mechanism will never be complete without the orbital physics taken into 
account. 
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